Two-particle scattering matrix of two interacting mesoscopic conductors 
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We consider two quantum coherent conductors interacting weakly via long range Coulomb forces. 
We describe the interaction in terms of two-particle collisions described by a two-particle scattering 
matrix. As an example we determine the transmission probability and correlations in a two-particle 
scattering experiment and find that the results can be expressed in terms of the density-of-states 
matrices of the non-interacting scatterers. 
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Recently there has been a growing interest in intrigu- 
ing subjects such as quantum measurement, controlled 
dephasing and shot noise correlations which involve two 
or more separate mesoscopic conductors interacting via 
long range Coulomb forces. We highlight here only three 
recent experiments. The noise cross-correlation of two ca- 
pacitively coupled quantum dots has been measured [l| 
and the sign of this correlation was found to be tunable 
by the gate voltage. In another work a Mach-Zehnder in- 
terferometer @, 0, S| coupled to a detector channel with 
shot noise was investigated both experimentally and the- 
oretically d, @] . Ref . Q shows how the current through a 
quantum point contact, capacitively coupled to a quan- 
tum dot, can be used to determine the statistics of charge 
transferred through the quantum dot. 

Several approaches have been developed to treat 
Coulomb coupled conductors. For systems in the tun- 
neling limit, a particle number resolved master equation 
approach 0, Q is often used. For systems well connected 
to contacts a self-consistent scattering approach has been 
developed and has been applied to dynamic conductance 
and charge fluctuations [l(J U, [H| in good agreement 




with experiment [13l | . A momentum resolved treatment 
of long range Coulomb interaction is necessary to treat 
the Coulomb drag which one conductor exerts on another 
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We are interested in a different formulation based on 
two-particle collision processes. We view two-particle in- 
teractions as the elementary process and derive a two- 
particle scattering matrix. In the weak interaction limit 
such an approach can be expected to have a wide range 
of applicability similar to Boltzmann equations with two- 
particle collision kernels. In contrast, in the strong in- 
teraction limit, it will be necessary to go beyond two- 
particle processes and permit for instance the interaction 
of a carrier in one conductor with a number of carriers 
in the other conductor @. 

Two-particle processes occur also if both particles are 
in the same conductor. Even in the absence of inter- 
actions, shot noise tests two-particle correlations and a 
predicted two-particle Aharonov-Bohm effect [l6( has re- 
cently been measured [171 ]. With interaction a number 
of highly interesting effects have been discussed in sys- 



FIG. 1: Two quantum dots coupled via the interaction 
\Q l Q 11 . The operators a' (c') create incoming electrons in 
the scattering states in the left (L) and right (R) leads of dot I 
(II), while b< (cP) are similar operators for outgoing electrons. 



terns with disorder [181 ]. for conductance and pumping 
fli^ ] and noise in quantum dots. In distinction, we 
extend scattering theory and obtain results that are both 
very general yet have still an immediate physical appeal. 

Interestingly the description of two-particle processes 
presented here involves, like the self-consistent scatter- 
ing approach , a gen eralization of the Wigner-Smith delay 
time matrix [2 ll. 1 2 21] . The Wigner-Smith matrix contains 
energy derivatives of the scattering matrix: its diagonal 
elements are proportional to the density of states gener- 
ated by scattering states incident from a particular con- 
tact [lOj. Thus diagonal elements are useful to describe 
the piled up charges in a conductor in response to a volt- 
age change at a contact. The off-diagonal elements of the 
matrix describe spontaneous charge fluctuations 1^, 11 1 . 

In Fig. Q] we show the type of system we are interested 
in. We consider two scatterers, both coupled to two non- 
interacting leads i = L,R. For simplicity we consider sin- 
gle channel leads. In the figure the scatterers are chaotic 
quantum dots (and we will refer to them as such), but 
our theory is valid for any scatterer. The Hamiltonians 
of the two (non-interacting) sub-systems read 



H + V 1 



H 11 = ff? 



H 



V 1 



(1) 



The symbols I /II refer to the first/second conductor. The 
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lead and dot Hamiltonians are given by 

Hj = E/ dEEa\(E)a t (E), H^Y^fifiU, (2) 

■i i 

Hj 1 = W dEEc\(E) Ci (E), H l d 1 = Y,9lg t E l .(3) 

i J i 

While a^E 1 ) creates incoming carriers in lead i of the 
first dot, c\(E) is a similar operator for the second dot. 
Furthermore fi (<?i) annihilates an electron with energy 
ti (Ei) in the first (second) dot. The coupling between 
leads and dots is described by 

yl = E / dE(4{E)fiWl+H.c), (4) 

il 

V U = J^f dE{c\{E) gi W% + H.c). (5) 

il 

We will assume an interaction of the form 

H c = ^Q l Q 11 . (6) 

with A a coupling energy and Q l = ej^ififi {Q n = 
e 12i 9i9i) tne charge operator on dot I (II). 

Before proceeding to treat this interacting problem, 
we recall properties of a single non-interacting dot (here 
dot I). The single-particle scattering matrix S I (E) relates 
operators which annihilate carriers in incoming and out- 
going states, hi = J2j^lj a j- The charge fluctuations at 
frequency lo = E — E' can be expressed in terms of a 
density of states matrix 

0,0 

V(E,g)-S*(K) g^jffi . (7) 

The diagonal elements of this matrix M\ L (and M^p.) 
are the injectance of the left (right) contact, i.e. the 
part of the total density of states associated with carriers 
incident in the left (right) lead. In the limit E — > E' the 
density of states matrix reduces to the famous Wigner- 
Smith delay time matrix discussed in the introduction. 

We are interested in correlations generated by the 
Coulomb interaction between the two dots. Consider the 
particles leaving the two dots. Particles in the outgoing 
states are created by operators b\(E), d\{E). The out- 
going two-particle state b\{Ei)d^{E 2 ) depends on two- 
particle input states in all pairs of incident channels and 
depends on the energy of the incident carriers not only 
at energy E\ and E 2 but also at E\ — u> and E2 + u> where 
uj is the energy which can be exchanged in the collision 
process. As a consequence the relation between incoming 
and outgoing two-particle states is 

bi(Ei)dj(E 2 ) = E [S I ik (E 1 )Sf l (E 2 )a k (E 1 )c l (E 2 )+ 

kl 

J duiSS ik ji(Ei,E 2 ,Ei -u,E 2 + u) 
a k {Ex-u))ci{E 2 +Lo)]. (8) 



Eq. ([8]) defines the two-particle scattering matrix. The 
first part of the relation is simply a product of the two 
single-particle scattering matrices, the second part con- 
tains the effect of the interaction. In this paper we will 
derive 8S up to first order in the interaction energy A. 

Let us first illustrate the derivation of the single dot 
scattering matrix (we choose dot I). We denote the n'th 
eigenstate of the Hamiltonian H d + Hj with energy E by 
(j) l nE . The eigenstates are either localized in the dot or 
lead. Since the spectrum in the lead is continuous, matrix 
multiplication also involves an integral over the energy. 
We are interested in the eigenstate %p l nE of H 1 which ap- 
proaches (fliE m the limit V 1 — > 0. It is expressed by the 
Lippmann-Schwinger equation [23| 

\^) = \€ E )+G I E k V I \^ E ), (9) 

with G l E = (E — H l d — H] ±ir)) . Here r\ is a posi- 
tive infinitesimal. It can be shown (see Ref. [23|) that 
I^Iie) (l^nis)) satisfies the boundary condition for ingo- 
ing (outgoing) states. The scattering matrix relates the 
two and is therefore defined to be Sj lm (E)S(E — E') — 
{i^nE^mE') ■ Using Eq. © one can rewrite it as 

SiJE) = S nm ~ 2mTl n (E, E), (10) 

with transition matrix T* m (E,E') = {<P l nE \V l \i^ EI ) [13]. 
For later use we will not only be interested in the tran- 
sition matrix between different lead states [24]] , but also 
between lead (1) and dot (d) states. Using the Lippmann- 
Schwinger equation we find 

T^ U \E,E') = W l [D\E')Y 1 W l \ (11) 
T^ ld \E,E') = W l [D\E')Y 1 (E' - Hj) , (12) 
T Kdl)( E ,E') = {E' - H l d )[D\E')Y 1 W l \ (13) 

with D l (E) = E-H\ + iirW^W 1 . Combining Eqs. (TTO]) 
and (jlip gives the well-known scattering matrix relating 
in- and outgoing lead states [IH 

S\E) = I 1 -2mW l [D l {E)Y 1 W 1 ^. (14) 

We will now turn to the case of two quantum dots. 
Without interaction (A = 0) the eigenstates are product 
states |^ Ei )®|V"e 2 ) = l<AU<£ 2 >- Let IC^,™^) be 
the eigenstate of the Hamiltonian H l ® I 11 + t l ® H n + 
H c which approaches V'mE,) if A — > 0. It fulfills a 
modified Lippmann-Schwinger equation 

\Cn El ,mB 2 ) = I ^mBi ^mE 2 ) + +E 2 He I Cn El ,mE 2 ) ' C^) 

with G E = {E~H l ® t n — I 1 (g) H u ±irf)- 1 . The two- 
particle scattering matrix has matrix elements 

S n mM{Ei, E 2 , E 3 , E i )5(Ei + E 2 — E 3 — E4) = 

(C nEl ,kE 2 \^mE 3 ,lE 4 )- 

(16) 
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With the help of the Lippmann-Schwinger equation, we 
find the interacting part 

SS nm .ki{Ei, E-2, E 3 , E4J — 

-2^l%^ Ei \H c \C nE3 , E4 ) = 
-2m(^g^ 1 \H e \^+ Es ^ l +) + 0(A 2 ).(17) 

For the last equality we have used Eq. (fT5|) and expanded 
up to first order in the coupling energy A. Because the 
coupling term H c is a direct product of operators working 
on dot I and dot II, we can write 



A 



We have expressed everything in single dot quantities and 
we can use the single dot Lippmann-Schwinger equation 
to proceed. The operator Q l is the charge operator of 
the dot and therefore 



1+ 



k 



E, 



(19) 



We only sum over states in the dot as indicated by the 
superscript (d). Using Eq. ([9|) we find 



{4e\^eJ = ^IeM 



TL(E,E 3 ) 



E% — E - 



irj 



(20) 



Combining Eqs. (T7]), (H]), © and ^ we find that 
the two-particle scattering matrix between different lead 
states depends on the single-particle T-matrices between 
lead and dot. Using Eqs. (fT2"|) and (fT3|) we calculate 



(21) 



5S(E 1 ,E 2 ,E 3 ,E 4 ) = ~2~ni\ 



gi(gi) - SHE 3 ) 
2ni{E 1 - E 3 ) 
S n (E 2 ) - S U (E 4 ) 
2iri(E 2 - Ei) 



Thus we have expressed the two-particle scattering ma- 
trix in terms of the scattering matrices of the uncoupled 
dots. Eq. (f2Tj) is the key result of this work. Its form 
reminds us of the density of states matrix Eq. ([7]). 

The application of our two-particle scattering matrix 
to a transport experiment such as the experiment in Ref. 
[l| requires that we take into account Pauli blocking in 
the Fermi sea of the leads. This is not trivial and we 
will present our solution in a later work 26]. Here we 
will illustrate the properties of the two-particle scatter- 
ing matrix by considering a real two-particle scattering 
experiment. We suppose our scattering problems to be 
one-dimensional and denote the coordinates along the left 
lead of the first (second) dot by x (y). We will assume 
that at time to a wave packet is created in each left lead, 
at positions xq and yo- The dots are at positions x = 



and y = 0. This means that we can describe the initial 
state in our system by 



I*) = J dEdE'a l (E)a ll {E')e 



i(k(E)x Q +k(E')y ) 



e- i( - E+E ^ to ^a{{E)c[{E')\0). (22) 

The functions a l {E) obey / dE\a l (E)\ 2 = 1, and k(E) > 
is the wave vector corresponding to an energy E. 

Since we work with wave packets the timing is impor- 
tant: if both wave packets reach the dots at very different 
times they cannot interact. Under the assumption that 
the width of the wave packet l/5k is much larger than 
vpTd, the dot appears effectively point like and we can 
factorize the influence of the interaction into a contri- 
bution from the wave packet overlap and a contribution 
from the scattering matrices of the dots. Here vp is the 
Fermi velocity and the dwell time in the dot. In this 
limit it is the integral 



2 = 



J dE 1 dE 2 dua I (E 1 )a I *{E 1 + uj)a n {E 2 )a n * (E 2 - lu) 



J(k(E 1 +Lj)-k(E 1 ))x +i(k(E 2 -uj)-k(E 2 ))y Q 



(23) 



that quantifies the overlap of the wave packets in 
the dot. For a Gaussian distribution a I / II (£') = 

e -(E-E F ) 2 /(46E 2 )/( 27T6E 2^1/4 w j th §E < Ep ( SQ ^ we 

can linearize the wave vectors around the Fermi wave 
vector hp) we find I = 2^/tt5E exp (— (x — yo) 2 &k 2 ). 
The widths in fc-space and -E-space are related by 5k — 
5E(m/2EFfi 2 ) 1 / 2 , with Ef the Fermi energy and m the 
mass. There is only an effect of the interaction if the wave 
packets are timed to reach the dots at about the same 
time, i.e. for \xq — yol^^ ^ 1 ( we assumed equal Fermi 
velocities). Furthermore X vanishes linearly with 8E, be- 
cause completely delocalized particles (plane waves) have 
a vanishing probability to be in the dot. 

Let us first calculate the probability that the first par- 
ticle leaves through the right lead, regardless of the be- 
haviour of the second particle. We define the operator 
n 1 R = j dEtf R (E)b[i(E) and we calculate 

(h l R ) = <*|nk|*>. (24) 



We rewrite the state (|2"2l in terms of output op- 
erators and d\ using the inverse of Eq. ([8]). 
Eq. ([24]) now contains expectation values of the 
form (01^(^)6,(^)6^(^)^(^)^(^3)4(^4)10) = 
S m S(E 1 - E 4 )6 jR S(E 2 - E)S Rm S(E - E 3 ). In the limit 
8k » vpTd (discussed above Eq. (f2"3")l ) the energy depen- 
dence of the scattering matrices can be neglected and 
we can evaluate all scattering matrices at the Fermi en- 
ergy. We will therefore suppress the energy argument. 
We write 



= T 



1 ^2( S RL S pL SS RL,pL 



T - IX 



dT 1 
~dE 



LL ■ 



Srl SpL $Srl ,pL ) 



(25) 
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FIG. 2: Quantum dot (system I) coupled to an edge state 
(system II). The quantum point contact (with transmission 
probability T and reflection probability R) is outside the in- 
teracting region. 



= IS^I 2 is the transmission probability of the 



Here T l 

first scatterer. Thus due to the interaction the transmis- 
sion probability of the first dot depends on the injectance 
of the left lead of the second dot defined in Eq. 0. 

Secondly we calculate the cross-correlation between 
particles in the right leads 



XI 



~dE 



LL 



' dE 



LL JV LR°LR) 



S^ErS^r) 



(26) 



with 5h 



i/ii 

R 



-I/II 



■I/H\ 



The cross-correlation de- 



pends on the off-diagonal elements of the density of states 
matrix of Eq. ([7]). For two completely symmetric scat- 
terers with Slr. = Srl and Srr = Sll it disappears. 

Let us work out Eq. (|26[) for two experimentally rel- 
evant systems. In a first example [H. 



27, 2S 



'J we 

assume that our second system is an edge state which is 
noisy because of the presence of a quantum point con- 
tact (QPC) left of it outside the interacting region, as 
shown in Fig. [2] (this is equivalent to assuming that the 
QPC has energy-independent transmission and reflection 
probabilities T and R = 1 — T). We do not specify the 
first system. For an interacting region of length L the 
density of states is TV 11 = mL/2irh 2 kp. We find 

(6h R 6h%) = -\I — N u RT. (27) 

The factor RT is a consequence of quantum partition of 
carriers at the QPC. 

With the experiment of Ref. [l| in mind, we now as- 
sume that system II is a quantum dot with a Breit- 
Wigner resonance at energy Ep = E n (so dT u /dE = 
at Ep) and with rates 7" (7") through the left (right) 
barriers. Conductor I is arbitrary. This gives 



7g) 



T(7i I + 7g)^ 



(28) 



Similar to the experiment the correlations Eqs. 
(126H28[) can have different signs. The sign of the correla- 
tion depends on the asymmetry of the scattering matrix 
but also on the sign of the energy derivative of the trans- 
mission probability [221 ]. 

To conclude we have calculated the two-particle scat- 
tering matrix for two weakly coupled mesoscopic con- 
ductors and we have expressed it as a function of the 
non-interacting scattering matrices. We illustrated the 
properties of this matrix by calculating the transmission 
probability and cross-correlation for a two-particle scat- 
tering experiment. Our results can be expressed in terms 
of the density of states matrices of the uncoupled conduc- 
tors. The approach developed here is very general and 
permits to treat a large class of systems. 

We thank M. Polianski and D. Sanchez for their valu- 
able comments. The work was supported by the Swiss 
NSF and the EU Marie Curie RTN "Fundamentals of 
Nanoelectronics" , MCRTN-CT-2003-504574. 
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